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Abstract 

The problem of magnetic vortex dynamics in an anisotropic spin liquid 
model is considered. For incompressible flow the model admits reduction 
to saturating Bogomolny inequality analytic projections of spin variables, 
subject the linear holomorphic Schrodinger equation. It allows us to con- 
struct N vortex configurations in terms of the complex Hermite polyno- 
mials. Using complex Galilean boost transformations, the interaction of 
the vortices and the vortex chain lattices (vortex crystals) are studied. By 
the complexified Cole-Hopf transformation, integrable N vortex dynamics 
is described by the holomorphic Burgers equation. Mapping of the point 
vortex problem to N-particle problem, the complexified Calogero-Moser 
system, showing its integrability and the Ifamiltonian structure, is given. 

1 Introduction 

As is well known the main ideas of Descartes vortex theory are represented in 
"Discours de la methode" (1637) and in a capital work "Principia Philosophiae" 
(1644). The Cartesian cosmology is based on primordial chaos, which by mo- 
tion according to fixed laws is ordering to cosmos [T]. According to Descartes 
the Universe is filled by thin all-penetrable fluid (similar to the ether) which 
is in a permanent rotational motion. The term vortex itself, " tourbillon" , is 
coming from comparison with turbulent motion of a river. Besides Descartes, 
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the vortex models of gravity were proposed by Bernoulli, Stokes and Huygens. 
But soon they were displaced by the Newton's gravity theory for a long time. 
Only in the middle of XIX century the interest to vortex theory revives with 
works of Helmholtz (1821-1894) [J, Thompson (Lord Kelvin) [3], and Kirchhoff 
(1824-1887) [1| on the vortex motion of an ideal fluid. The mathematical de- 
scription of processes related with the motion of vortex in a liquid is starting 
from Helmholtz's paper "Uber Integrale der Hydrodynamischen Gleichungen 
Welche den Wirbelbewegungen Entsprechen" (1858) in which he formulated his 
famous theorem on conservation of vorticity in the rotational motion of a fluid. 
He also notices an analogy between the fluid motion and the magnetic action 
of electric fields. General equations of motion for N point vortices (Kirchhoff 's 
equations) have been introduced by Kirchhoff in his lectures in mathematical 
physics |4]. He derived the corresponding Hamiltonian form of equations and 
found all possible integrals of the motion. In contrast with Newton's equation 
for A'" point masses, having the second order, the Kirchhoff equations are the 
first order of vortex coordinates. Walter Grobli in his thesis "Specille Probleme 
uber die Bewegung Geradliniger Parallerer Wirbelfaden" in 1877 [5j analyzed 
the integrable problem of motion of three vortices in the plane. He obtained 
the system of three nonlinear equations possessing two integrals of motion and 
allowing to get explicit quadrature. He also considered particular case of the 
problem of four vortices under condition of symmetry axis and more general 
problem of 2N vortices with N symmetry axes. 

The interest to vortex theory increases with Kelvin's vortex theory of atoms 
"On Vortex Atoms" , [3j. But soon his model was dismissed by quantum me- 
chanical model of atoms [6]. Kelvin also posed the problem of stability under 
the stationary rotation of the system of N point vortices located at N polygon 
vertices. He noticed that the problem is similar to the problem of stability for 
the system of equal magnets floating in external magnetic field. Experiments 
with floating magnets performed by Mayer, leaded him to the conclusion that 
for number of vortices (magnets) exceeding 5, the rotating polygon becomes 
unstable (in fact the case iV = 6 is stable). The linear stability of polygon has 
been studied by J.J. Thompson. He found that for n < Q the linear stability 
takes place, while for rt ^ 8 is not. Stability of the case n — 7 needs nonlinear 
analysis and after several failed attempts has been proved only recently. Non- 
integrability of the four vortex problem in plane, indicating on chaotization of 
the vortex motion, was found recently [7]. 

Many problems involving interfacial motion [8] can be cast in the form of 
vortex sheet dynamics ^9j. The discovery of coherent structures in turbulence, 
increases expectations that the study of vortices will lead to models and an 
understanding of turbulent flow, one of the great unsolved problems of classical 
physics. Vortex dynamics is a natural paradigm for the fleld of chaotic motion 
and modern dynamical system theory [TU] . The theory of line vortices and vor- 
tex rings is a part of modern theory of liquid Helium II. Interaction of vortex 
structures essentially influences on processes in atmosphere and the ocean. In 
techniques, complete understanding of friction to the motion, noise generation 
and shock waves, is impossible without clear theory of the vortex motion. In 
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a wide and important class of motions of ideal inviscid fluid [TT], the vortex 
dynamics provides physically profound examples of nonlinear Hamiltonian sys- 
tems of infinite dimensions, attracting much interest in relation with chaotic 
phenomena in dynamical systems [12] . 

Modern applications of vortices extend from liquid crystals and ferromagnets 
[IB] , [Hj to superfluids [T^ , [H] and from non-equilibrium patterns to Quantum 
Hall effect [17] and cosmic strings [18] , [19] , [20] . They can play essential role in 
the loop gravity and quantum spacetime, leading to formation of the Cantorian 
Space [3T]- fractalization of the microspace geometry, which has been discussed 
recently in the context of the fractional Quantum Hall Effect [55], [23] and 
Cosmic Strings [21] ■ 

Recently it was shown that in ferromagnetic nanomagnets, for particles of 
various shapes (strips and rings), the switching process involves domain walls 
which are composite objects made of two or more vortices and edge defects 
with integer or fractional winding numbers [25j . It may provide a basic model 
of complex magnetization dynamics in nanomagnets by reducing it to the cre- 
ation, propagation, and annihilation of a few planar topological defects [5S] . 
The planar vortices represent important class of magnetic systems admitting 
topologically nontrivial solutions, the dynamics of which is one of the most in- 
triguing questions. Due to nonlinearity in general it is very difficult task to 
construct exact analytical solution even for one vortex. This is why simplified 
models admitting exact treatment are important tools to study more realistic 
situation. In the present paper we consider the simple model of anisotropic pla- 
nar spin liquid, admitting construction of exact N-vortex and N- vortex lattice 
solutions with integrable dynamics for arbitrary N. 

2 Topological Spin Liquid Model 

Inspired by superfluid motion at T = in He^ — A, the model introduces the 
normal velocity of the fermionic liquid as an additional hydrodynamical variable, 
describing the background fermionic vacuum. The magnetic fluid or spin liquid 
is characterized by the local magnetization field S{x,y,t) subject to the modi- 
fied Heisenberg model in the moving frame with velocity v{x, y, t). Moreover, for 
planar magnetic systems the existence of topologically nontrivial vortex config- 
urations requires the fluid to be rotational with non- vanishing vorticity function 
[55] . The system is 

St + vidiS - v2d2S = S X {df - di) ■ S (1) 

diV2 - d2Vi = 2S{diS xd2-S) (2) 

where S'^{x, y,t) — 1 is classical spin field, v — (wi, ^2) is velocity field. The first 
equation is anisotropic in x and y directions, the second equation, requiring that 
the fiuid vorticity in the plane is proportional to the corresponding magnetic 
topological current component, is known in the theory of superfiuid He^ as 
the Mermin-Ho relation |15j . [16] . This system in general should be supplied 
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with the continuity equation pt + W ■ {pv) — for the density p{x,y,t). For 
incompressible flow pt = the last equation simplifies to V • t/ = and allows 
one exclude p from consideration. Moreover, for the fluid flow constrained by 
the incompressibility condition 

divi + d2V2 = 0, (3) 

the conservation law 

dtJo + d2J2-diJi = (4) 

holds [35], where 

Jo = {dj)^ + {d2S)\ (5) 
Ji = -2diS-Sx{d^-di)S + viJo + 2S-{diSxdiS-did2Sxd2S) 
J2 = 2d2S-Sx{d^-d^)S + v2Jo-2S-id'lSxdid2S-diSxd2S). 

Due to this, for the incompressible flow ^ the energy functional 

E= f [ Jod^x = [ [ {{diSf + {d2Sf}d^x (6) 



is the conserved quantity. The topological charge or the winding number of a 
spin configuration is defined as 



1 
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S ■ {diS X d2S)d^x (7) 



and it is also the conserved quantity. These two quantities are related by the 
Bogomolny type inequality 

E > 87r|Q|. (8) 

It follows from the evident one 

{d^S±eijS X djSfd^x > (9) 

and is saturated by time dependent spin configurations satisfying the self-duality 
equations 

d,S±eijS X djS ^0 (10) 

If the spin vector phase space, the 2D sphere, we consider as the Riemann sphere 
for a complex plane, we can project points on this sphere to that plane by the 
stereographic projections 

2C 1 - ICP 

s+ = Si + iS2 = ^ ^ |^|2, s'3 = Y+]cF' 

where (^{x,y,t) is a complex valued function. Now we rewrite the self-duality 
equations PH)) in stereographic projection form. By the complex derivatives, for 
the first sign in (|10[) we have the analyticity or the holomorphicity condition: 

C.-(a^,y,i) =0 (12) 
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while for the second sign we have the anti-analyticity or the anti-holomorphicity 
condition: 

C.(x,y,i)=0 (13) 

where z — x + iy and dz — ^(d^ — idy), dg ~ ^{d^ + idy). These conditions 
written in terms of the real 3?^ and imaginary 3^ parts of function C, represent- 
ing the Cauchy-Riemann equations, describe the incompressible and irrotational 
fluid flow with the velocity potential and the stream function 3^ [27] , 



3 Anti-Holomorphic Reduction 

As we can see, analytic/anti-analytic configurations saturate Bogomolny in- 
equality and have minimal energy. This suggests to find solutions of the model 
([T]) with holomorphic/anti-holomorphic stereographic projections. For this rea- 
son we first rewrite ([1]) and Q in stereographic projection form in complex 
coordinates z for complex velocities v± = vi ± iv2 and if C is anti-holomorphic 
Cz = 0, then the system is reduced to 

iCt + iv+C, + - 4^^^ = 0. (14) 

dzv+ - d-zv^ = -Uj^^^^. (15) 

To be consistent, this anti-holomorphicity constraint must be compatible with 
time evolution. By direct computation we can show that for incompressible flow 
([3]) the anti-holomorphic constraint = is compatible with the time evolution 
dC,z/dt^{). For this flow from ^ and US]) we have Ct + FCs = and = 0, 
where F = vj^ — 2i In . 



4 Anti-Holomorphic Ishimori Model 

The above consideration suggests to solve the incompressibility conditions ex- 
plicitly. So we consider the topological magnet model ([1]) , ([2]) with incompress- 
ibility condition Equation \7 ■ v — can be solved in terms of the stream 
function of the flow, vi = 8211^, V2 — — so that we get the so called Ishimori 
Model [ini- 

St + 82^818 + di^pd2S = Sx {8fS-diS) (16) 

{df + 81)11; = ~2S- {818x828). (17) 

The Ishimori model is the first example of integrable classical spin model in 
2-1-1 dimensions [3S]. It was shown to be gauge equivalent to the Davey- 
Stewartson equation, representing the 2-|-l dimensional generalization of the 
Nonlinear Schrodinger equation [32], [S3], [37] • Though it was solved in terms 
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of the d problem, for description of vortices and vortex lattices we propose more 
simple and elegant method. In terms of complex variables 

v+ = vi+ iv2 = -2iipz, V- = vi- iv2 = 2i^pz (18) 

for incompressible flow, preserving anti holomorphicity constraint, we have de- 
pendence C = C{z, t) and the model reduces to the system 

+ 2V',C.- + - Y^i^ = (19) 

We can rearrange the first equation as 

<t + Kd^ - 21n(l + + InC.:}.- = 0. (21) 
If we choose = 2 ln(l + then equation pp)) 



1 + ICP 



is satisfied automatically. Then from (PT|) for function we have the anti- 
holomorphic time dependent Schrodinger equation 

<t + - 0. (23) 

By complex analog of the Cole-Hopf transformation u = 4:(z/C = 4(LogC)j it 
implies complex Burgers' equation 

iut + uuz + 2uzz — 0. (24) 

In u can be interpreted as the complex velocity of effective flow with the 
complex potential f{z) = LogC,'^. Then every zero of function C, corresponds 
to the vortex solution (pole of complex velocity) of anti-holomorphic Burgers' 
equation. 



5 Magnetic Vortices as Moving Zeroes of Her- 
mite Polynomials 

By stereographic projection pT|) at every zero of function ({zk,t) = we have 
iSi+iS2)izk,t) = 0, 53(zfc,i) = l (25) 
From another side for the polynomial Cjv of degree N at infinity |z| — s- cx) 

iSi+iS2)izk,t)=0, S3izk,t)^-1 (26) 
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Figure 1: iV = 1 Magnetic Vortex 



It shows that every zero corresponds to the magnetic vortex located at that 
zero with the spin vector S directed up, while at infinity it is directed down 
(ferromagnetic type order)(See Fig. 1). With such boundary conditions the 
topological charge Q in ([7]) is an integer valued and characterizes the number 
of magnetic vortices. 

If we calculate the topological charge ([7]) for N zeroes solution 

N 

C{z,t) = l[{z-zk{t)) (27) 

fc=i 

then we find that Q — —N. The above consideration shows that motion of 
zeroes of time dependent Schrodinger equation (|^ describes the motion of 
magnetic vortices in the Ishimori Model. Moreover due to (|24p dynamics of 
these vortices is described by the anti-holomorphic Burgers' equation. To study 
motion of zeroes of ([^5]) we construct the generating function of the basic vortex 
solutions by considering complex plane wave solution 

C(z,t) = e'=^+2''='*. (28) 

Let X = k^J^, then we rewrite it as the generating function for the Hermite 
polynomials of complex argument 



or 



gfc^+2^fc^^ ^ g-x^+2(.V^)- = V U J - I ^ (29) 



n— 




C(z,t) = ^ -{-2ttr/'H,, ( zJ± ] = J2 ^*»(^"'*) (30) 



where at every power k" we have a polynomial solution of order n: 




(31) 
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This polynomial has n complex roots zi{t), z„(i) describing positions of vor- 
tices. For complex zeroes of this function, by identification 

^Niz,t) = l[{z-Zk{t)) = {-mr/^H^ ( ^^^^ 

we find that z = Zk{t) implies Hn{j^^^) = 0. Denoting w^^\ k = 1, ...,n as 

zeroes of Hermite polynomials, Hn{w^^) = , we have the time dependence 
for vortex positions 

Zk{t) = 2w[^^yf^. (33) 

Due to reality of roots wlf ^ the form of our solution Zk{t) = 2w'^^'' \/ —2it 
implies that all vortices are located on diagonal lines of complex plane: 

Zk{t) = |2«;f ^v^le^"/^ (34) 

Wc note that since the time dependence includes square root of time variable 
t, then under time reflection, when t is replaced by —t, position of vortices will 
rotate Zk — > e^f^Zk on angle 7r/2. It means that under collision our vortices 
change velocity in orthogonal direction and from one diagonal line would be 
displaced to the orthogonal one. Moreover, sum of vortex positions, representing 
the center of mass of the system, is integral of motion located at the beginning 
of coordinates. 
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In Figure 2 we show contour plot of basic vortex dynamics for four vortices. 
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Figure 2: N = A Vortex Dynamics 



In Figure 3 interaction of two magnetic vortices in Ishimori model is shown. 




Figure 3: N = 2 Magnetic Vortex Dynamics 



Since function ^f^v in is solution of the linear equation for any 
integer N, then any linear combination of these functions 

$Ar(z, t) = aAr*Ar(z, t) + aN-l'i' N-l{z, t) + ... + ao*o(^, t) (35) 

aN ^ is also a solution. This solution is determined by N complex constants 
ai, 02, ajv, which are integrals of motion of the system. The higher order 
coefficient oat 7^ is not essential and could be put to one. Below we represent 
particular cases for A'' = 3 and N = A. 

1) For N = 3 we have the general solution 



t) = (z^ + 12zit) + 02(2^ + -iit) + aiz + cq 



(36) 
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This cubic in z equation has 3 complex roots z\{t\ Z2{t), z^it) moving in plane 
according to the systems (|76l) . Instead of solving that system of differential 
equations we will find roots of cubic equation according Cardano formulas. Co- 
efficient 

-aa = zi{t) + Z2it) + zsit) 

is integral of motion having meaning of the center of mass for three vortices. 
Without loss of generality we can always choose coordinate system with the 
beginning at this center of mass. So we will put 02 = 0. Then our cubic 
equation has the reduced Cardano form 



z 
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pz + q^O (37) 



where p = ai + \2zit, q ^ ao + 4a2«i. Solution of this equation is 

zi = ai+/3i, Z2 ^ aiUJi + P1UJ2, Z3 ^ aiUJ2 + PiUJi (38) 
where ai,/3i is one of the couple of roots 

/ I \ 1/3 

(39) 




and 



ivi = -1/2 + iV3/2 = e^2V3^ ^2 = -1/2 - iV3/2 = e'^^Vs 



are cubic roots of 1. For particular values ao — a2 — when q = our roots 
become 

= Vliiil + e'"/^) (40) 
Z2 = V^t{uji+uj2e'''/^) (41) 
Z3 = V4^^(L02 + cjie"/3) (42) 



and coincide with particular cases (I33p when one of the vortices, Z2 = 0, is static 
at the beginning of coordinates. In general case due to non-vanishing orbital 
momentum no one of three vortices crosses beginning of coordinates. 
2) For N = 4 case we have 

$(£,<) = {z^ + 2Az'^it+\2{~2itf) + a^{z''^ + \2zit) + a2{z'^+Ait)+aiz + ao (43) 

and like in previous case we choose the center of mass at the beginning of 
coordinate system, so that the coefficient 

-03 = zi{t) + Mt) + Mt) + zi{t) = 

Explicit form of the roots in this case can be done in radicals. Below we will 
write them only for special case ai = 0. Solving bi-quadratic equation z'* -|- 
z^[a2 + 2Ait] -I- [ao -|- 4a2ii -f 12(-2it)^] we have four roots: 



21,2 = ±1 



^ + mt 



y + 12ii) - [ao + 'ia2it + l2{-2ity] (44) 
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Z3,i = ±\l- ^ + 12it - + 12ii)'- [ao + 4a2ii+12(-2it)2] (45) 

From dynamics of 3 and 4 vortices we can see that differences appear for the 
vortex motion at finite times. But for large time the behavior of vortices is 
similar to the particular case ([33]) . Explanation of this fact can be done for 
general N-vortex configuration. Indeed, if we consider asymptotic form of the 
general N-vortex solution ([SS]) . when t ^ oo and z — > oo such that \z\'^/t 
const, then we can see that all terms of the function for any N have the 
same order. Moreover, the dominant role in ([55]) plays the function '^n with 
highest order of N. But it is exactly solution ([5^ which we found before. So we 
proved that asymptotically our vortices will follow diagonal lines according to 
the low ([55)1 : Zk{t) — 2wkV—2ii. If we calculate complex velocity corresponding 
to k- vortex Uk = dz^/dt = Wk\/—2i/t then at large times t —f oo the velocity 
of the vortex is decreasing up to zero Ufe ^ 0, as the inverse square root of time 



6 Vortex Solutions Generating Technique 

In this section we propose a general method allowing us to create an arbi- 
trary number of vortices on given background solution. In general, N vortex 
configuration is described by complex polynomial function degree N where co- 
efficients are represented in terms of symmetric polynomials. The polynomial 
Pi = zi + .. + zn is integral of motion having meaning of the center of mass 
for TV vortices. However, all other polynomials are not integrals of motion, only 
proper combinations of these polynomials provide integrals of motion of the sys- 
tem. Moreover, to find solutions zi{t), ...,2^(1) in terms of these polynomials 
means solving algebraic equation degree N, which as known to be solvable in 
radicals only for < 4. This is why to add a new zero or vortex to the system 
we follow another approach. In section 4 we have constructed basic solutions 
in terms of Hermite polynomials of complex argument (|32p . Adding recurrence 
relations 

HN+ii^) = 2xHNix) - 2NHn-i{x) (46) 
= -H'jfix) + 2NHn-i{x) (47) 

one can get 

Hn+i{x) = {2x - ■^)Hn{x) (48) 
dx 

Recursively it leads to following operator representation of the standard Hermite 
polynomials [38] 

HN{x) = i2x-^f -1 (49) 
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which implies for complex argument 

^Niz, t) = (-2»i)^/2(-i= - 2V^t^.f ■ 1 (50) 
and we have operator representation for our basic solution (j32p 

^N{z,t) = {z + m^)'' -i (51) 

az 

This representation suggests the form of generating operator for solutions of our 
equation (j23[) . By direct substitution we can prove the next, vortex generation 
technique for solutions of this equation. If <i>(z, i) is a solution of pS)) . then 
function 

^(z,t) = (z + Ait^)^z,t) (52) 

is also a solution of (221). This solution add one vortex to the background 
configuration t). Applying it several times one can add arbitrary number 
of vortices. Equation (I23p has evident solution $ = 1. Then 

^-1(2, t) = {z + Ait—)^{z, t) = {z + Ait—) ■ 1 (53) 



is also a solution of (|23|) . Next we have 

^2{z, t) = {z + 4^^— )*i - (z + 4zt— )2 . 1. (54) 
oz oz 

Continuing this procedure we have solution of (|23p for an arbitrary positive 
integer N 

d d 
^n{z, i) = (z + 4^^— )*Ar_i = ... = {z + 42t— )^ • 1. (55) 
Oz oz 



This way we derived particular operator representation for N vortex solution 

N Q 

*Jv(z, t) = Y[{z - z^{t)) = iz + ^it—fl. (56) 



Using previous result and linearity of we have the next generalization. If 
$(z, t) is a solution of equation ((23l) then function 

N , 

vl/(z,i) = ^a„(z + 4zi— )"<i>(z,t) (57) 

n=0 



where ao, ajv are arbitrary constants, is also a solution of equation (j23p . As 
easy to see adding to the system a new vortex in a proper way, we add an 
additional integral of motion. Let us suppose that we have solution '^{z, t) with 
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N - simple zeroes at points zi, ...,zn, that means, ^'(z„,i) = 0, n — 1, ...,7V. 
Then due to 

N N 

*(^, t)=Y[iz- zn{t)) = a„*„(z, t) (58) 

n— 1 n—0 

we have the system of N linear algebraic equations X]^=o '^n'^n{zk, t) = 0, fc = 
1, ...,N. Extracting n — term and dividing on gq it can be rewritten in the 
form of inhomogeneous system of N algebraic equations 

N 

J2bn'in{zk,t)^-1, k = l,...,N (59) 

n=l 

on N variables &„ = a„/ao. Then, N integrals of motion can be found by 
Crammers formulas bk = A^/A, fc = 1, ...,N. From the above consideration if 
we have <l>(z, t) as a solution of equation (|23p and -F'(z) is anti-analytic function 
in some domain Dq — {\z\ < i?}, then function 

*(z,i) = F(z + 4it-^)$(z,t) (60) 

is also a solution of equation ([23|) . Since coefficients a„ = ''^ ^^/"-^ are integrals of 
motion, the function F can be considered as the generating function of integrals 
of motion. At the end of this section we note that operator K = z + Aitdg in 
our approach is commuting with the Schrodinger operators' = idt + 29| and 
represent the complex boost transformation of the Galilean group. 



7 Sine-Hermite Solution and Stationary Vortex 
Lattice 

As an application let us consider entire function 

c=o 2n+l 

F(u.) = sinw; = y(-l)"y^^^ -. (61) 

^ ' ' 2n-|-l)! ^ ' 

n—0 ^ ' 

Then according to the above consideration 

n—0 ^ ^ ^ 

is solution of equation According to (fSTjl we rewrite it as 

n—0 
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or by ([5^ 



(-1)" 

(2n+ 1)! 
1 



-(-2zi)"+^iJ3 



1 



{-2it)"+2Hi 



+ -(-2*t)"+5i/5 
5! 



2V^2it 



Writing operator 



sin z + 4ii 



1 

2i 



exp 



i [ z + at 



dz 



exp 



i [ z + at 



dz 



(65) 



and using Baker-HausdorfF relation for non-commuting operators z, we have 
d 



or 



sin z + 4ii 



dz 



dz 



d 



d 



sin I z + 4it— I —ie * [cosz sinh(4i— ) — isinz cosh(4t— )] 



dz 



dz 



Applying this operator to evident solution '5 = 1 we get 



d 



sin 1 z + 4:it— ] ■ 1 = e smz. 
dz , 



(66) 



(67) 



(68) 



It describes the time dependent stationary vortex lattice solution of (^5]) . Po- 
sitions of vortices in the lattice do not change with time. Comparing with ((62] 
we have the next expansion in terms of Hermite polynomials 



i(-2ii)"+ii73 



exp(-2ii)sinz = (-2ii)"+2i7i 



1 



2V-2it 



,2V^2i</ 5! 
Then in terms of it gives 

$(z, t) = e~^** sin z = ^i (z, t) 



2y/-2it 



3! 



5! 



(69) 
(70) 

(71) 



The last formula describe decomposition of the stationary periodic lattice in 
terms of motion of odd l,3,5,...,2n-|-l,... number of vortices. It shows that 
proper superposition of odd number basic vortex motions described in Section 
5 leads to the fixed in time periodic lattice of vortices (vortex crystal). 
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8 Single Vortex - Vortex Lattice Collision 

As another application of the resuhs from Section 5, now we construct a new 
class of solutions describing single vortex collision with the vortex chain lattices. 
Consider solution of in the double lattice form 

C(z,t) = e"^"sin(z-zi(t))sin(z + zi(i)), cos 2zi (t) = re^** (72) 

then we have another solution of (|23p in the form 

d 1 

= (z + 4it— )-(r-e~***cos2z) 
oz 2 

= i[re^**z-zcos2z + 8iisin2z]e"^'* (73) 

Using properties of trigonometric function of complex argument 

sin 2z = sin 2x cosh 2y — i cos 2x sinh 2y (74) 

we have for the real and imaginary parts of function ^' following expressions 
correspondingly 

r r 11 

3?^' — -a;cos8i+-'i/sin8t+8t cos2a;sinh2-;/ a; cos 2x cosh 2t/ — w sin 2a; sinh 2-;; 

2 2 2 2 

r r 11 

= -X sin 8t— -y cos 8i+8t sin 2x cosh 2y+-y cos 2x cosh 2y— -x sin 2x sinh 2?/ 
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In figure 4 we show collision of a single vortex with the double vortex lattices 
at positive time i > 0. As we can see, addition of the vortex leads to dimerization 
of the lattice vortices which propagates in both directions and creates finally 
the dimerized lattice. 




Figure 4: Single Vortex - 2 Vortex Lattice Dynamics 
Following the same idea we can consider solution describing interaction of 
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N- vortices with M- vortex chain lattices (vortex crystal) in the form 



M 



^{z, t) = e-2'*^'*(z + Att^_)^ Y[ sm{z - Zk{t)). (75) 

9 Integrable N-particle Problem for N- Vortex 
Motion 

In this section we show that problem of N-point vortices in the plane can be 
reduced to complexified version of the Calogero-Moser model type I [30], [3T] . 
As we have seen in section 4 the system of N point vortices is described by 
function ([27]) satisfying the anti-holomorphic Schrodinger equation psp. Then 
positions of vortices in the complex plane, zi,...,zjv, are subject to the first 
order system 

d ^ I 

In one space dimension this system has been considered first in [32] j [30] for 
moving poles of Burgers' equation, determined by zeroes of the heat equation. 
However, complexification of the problem has several advantages. First of all 
the root problem of algebraic equation degree N is complete in the complex 
domain as well as the moving singularity analysis of differential equations. In 
contrast to the one dimension, the pole dynamics in the plane becomes time 
reversible (see discussion after eq. (IMt ) and has interpretation of the vortex 
dynamics. If we differentiate (|76p once and use the system again then we have 
the second order Newton's equations of motion 

These equations have Hamiltonian form 

dH . dH 

with Hamiltonian function 

J = l ]<k ■' 

The system (|771) is complexified version of the Calogero-Moser system type 
I, where N-particle positions, qi,...,q]\[ should be replaced by complex vortex 
positions Zi,...,Z]\[. The Hamiltonian equations (j78[) are equivalent to the Lax 
matrix equation 

iL = AL- LA (80) 
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where 



Ljk = SjkPj + i.g(l - Sjk)- 



Zk 



Ajk 



[Z] -zif 



(1 



(Zj - ZkY 



(81) 



(82) 



and the couphng constant g = \/— 4. Since matrix L{t) is isospectrally deformed 
with time, the corresponding (complex) eigenvalues are time independent inte- 
grals of motion. Their symmetric functions as integrals of motion, are given by 
Ik = trL^'^^ . It shows that complexified Calogero-Moser system is an integrable 
system and as a consequence, the N- vortex system (j76p . which has been mapped 
to Calogero-Moser system, is also integrable. 



10 Integrable N-particle Problem for N- Vortex 
Crystals 

Similarly to the previous case now we consider mapping of the N-vortex chain 
lattices (vortex crystals) to the complexified Calogero-Moser system of type II 
and III |31| . For simplicity first we consider the system of two vortex chain 
lattices described by function 

C{z, t) = e"'^'* sin(z - zi{t)) sm{z - Z2{t)) (83) 
so that position of lattices is subject to the first order system 

zi — — 4icot(zi — Z2), Z2 — 4icot(zi — Z2) (84) 
Differentiating once in time we get Newton's equations in the Hamiltonian form 



dH . dH cot zi - Z2) , , 

zi^ T— ^Pi, Pi = --^ = 32 . — (85) 

dpi dzi sm''(2:i - Z2) 



dH . dH cot(z2 - zi) 
Z2 = ^ = P2, P2 = = 32 . — (86) 

dp2 dz2 sm'^(z2-2;i) 

with Hamiltonian function 

z z sm [zi — Z2) 

Comparing this Hamiltonian of two vortex lattices with the Calogero-Moser 
system, we realize that it corresponds to complexified version of the model type 
HI. We can generalize this result considering N vortex chain lattices in the 
horizontal direction x. Positions of lattices are subject to the first order system 

N 

^3 = - 3 = 1, N. (88) 

fe=i,(fc#j) 
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Differentiating once we get 



I,. -32 £ n^-^^ .^l,...,iV (89) 



which is complexified Calogero-Moser system type III with Hamiltonian 



If instead of horizontal x direction, we consider N chain lattices in the vertical y 
direction, it results in rotation of every zero on angle 7r/2 and the replacement 
of complex function sinz by sinhz. As a result, the corresponding Calagero- 
Moser system would be of type II. This consideration shows also equivalence of 
complexified Calogero-Moser systems of type II and III. 

11 Time Dependent Schrodinger problem in 
Harmonic Potential 



Vorticity equation (|22p is invariant under substitution ^ tp + U where U is 
an arbitrary harmonic function: AU — 0. Choosing 

V' = 21n(l + |Cp) + t/(z,t) + C/(z,t) (91) 

and substituting to (|2ip we have complex Schrodinger equation with additional 
potential term 

Kt + C.-.- + C,S-z - (92) 
12 Bound State of Vortices 

Here we choose particular form i) = ^ so that 

^ = 21n(l + |C|') + i(z2 + z). (93) 
Then we have time evolution subject to the equation 

+ K-z-z + K-z - 0. (94) 

Looking for solution in the form 

C(z,i)=E^"*""(^) (95) 

n 

we find that function Uniz) satisfies the complex Hermitc equation 

u" + ZM^j + M„ = 0. (96) 
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It gives time dependent vortex solution in the form 

N 

C(z,t) = ^e™*i/„(z). (97) 

ri=0 

This solution is A'^-th degree polynomial with periodic time dependent coefR- 
cients and it describes the bound state of N magnetic vortices. For particular 
value N = 2 we have 

C{z,t) = Ho{z) + e''Hi{z) + e''''H2{z), (98) 

5RC = 1 + 2xcosi + 2ysint + [4(a;^ - j/^) - 2] cos 2i + Sxy sin 2t (99) 

3C = -2ycosi + 2xsini- 8xycos2i + [A{x^ - y^) - 2]sin2t (100) 

This solution is periodic in time with period T — 271 and it describes the bound 
state of two magnetic vortices (see Fig. 5). 




-2 -1.5 -1 -0.5 0.5 1 1.5 



Figure 5: Bound State of Two Magnetic Vortices 



13 Static Vortex Crystals and the Liouville Equa- 
tion 

From previous considerations we have seen that static vortex configurations for 
Eq. ((23|) are exist only for N = 1. One can ask if there exist static N vortex 
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or even N- vortex lattices- vortex crystals for the Ishimori model (|16p . p?)) . To 
answer this question we return back to ((2T|) . For static configurations d(/dt = 
implies 

^ = 2 ln(l + - In G - In C. = - In J^T^ ^^^^^ 

so that the vorticity equation ((22)) is satisfied automatically. But we notice that 
pOip is the general solution of the Liouville equation 

AiP = Se"'^. (102) 

If one chooses C = Ili^il^" ^0 with N simple zeroes in the complex plane then 
it determines N static vortices located in the plane at zeroes of this function 
[34j . To have static vortex lattices periodic in x direction we can consider 

N 

C^Y[sm{z-z,) (103) 

It determines N time independent magnetic vortex lattices. For = 1 the 
stream function 

+ ^ (104) 

I coszj"' 

describes periodic in x static lattice of vortices. If we choose ^ = zsinz then 
the stream function 

I sm z + z cos 
describes single vortex imposed on vortex lattice. 



14 CONCLUSIONS 

To solve the problem of magnetic vortices in a planar spin liquid model, first 
we found holomorphic reductions of the model and showed that evolution equa- 
tion at this reduction becomes the linear anti-holomorphic Schrodinger equa- 
tion. Analogy of the anti-holomorphic Cole-Hopf transformation with the well 
known hydrodynamical relation between the complex velocity and the complex 
potential, leaded us to formulation of the complex Burgers' equation with inte- 
grable N vortex dynamics. We found that vortices correspond to zeroes of the 
complex Schrodinger equation. This allowed us to construct N vortex configu- 
rations, vortex chain lattices and their mutual dynamics in terms of the complex 
Hermite polynomials. By mapping our vortex problem to N-particle problem, 
the complexified Calogero-Moser system, we showed its intcgrability and the 
Hamiltonian structure. 

Finally we note that the holomorphic Hopf equation 

iut + uuz = (106) 
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which corresponds to the dispersionless limit of the holomorphic Burgers' equa- 
tion (|24p . has been considered very recently as nonlinear bosonisation in quan- 
tum hydrodynamics for description of quantum shock waves in edge states of 
Fractional Quantum Hall Effect [ID]. The weak solution of this equation for 
point vortices with strength Fi, ...jFat, so that 



When all the vortex strengths are equal Fi = ... = Fjv then this system reduces 
to ((76| and is integrable. However, in the general case the system is not known to 
be integrable. In particular, for = 3 the system with constraint Fi = F2 7^ F3 
has been studied in [41j to explain the transition from regular to irregular motion 
as travel on the Riemann surface. 
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